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Abstract — We consider the problem of decentralized power 
allocation for competitive rate-maximization in a frequency- 
selective Gaussian interference channel under bounded channel 
uncertainty. We formulate a distribution-free robust frame- 
work for the rate-maximization game. We present the robust- 
optimization equilibrium for this game and derive sufficient 
conditions for its existence and uniqueness. We show that an 
iterative waterfilling algorithm converges to this equilibrium un- 
der certain sufficient conditions. We analyse the social properties 
of the equilibrium under varying channel uncertainty bounds for 
the two-user case. We also observe an interesting phenomenon 
that the equilibrium moves towards a frequency-division multiple 
access solution for any set of channel coefficients under increasing 
channel uncertainty bounds. We further prove that increasing 
channel uncertainty can lead to a more efficient equilibrium, 
and hence, a better sum rate in certain two-user communication 
systems. Finally, we confirm, through simulations, this improve- 
ment in equilibrium efficiency is also observed in systems with 
a higher number of users. 

Index Terms — Game theory, rate-maximization, Nash equilib- 
rium, waterfilling, CSI uncertainty, robust games. 



I. Introduction 

Rate-maximization is an important signal processing prob- 
lem for power-constrained multi-user wireless systems. It 
involves solving a power control problem for mutually in- 
terfering users operating across multiple frequencies in a 
Gaussian interference channel. In modern wireless systems 
where users may enter or leave the system freely and make 
decisions independently, decentralized control approaches and 
distributed algorithms are necessary. Game-theoretic methods 
provide an appropriate set of tools for the design of such 
algorithms and have been increasingly used for the analysis 
and study of communications problems 0. 

In multi-user systems, the users can either cooperate with 
each other to achieve a socially optimal solution or compete 
against one another to optimize their own selfish objectives. 
Cooperative game-theoretic approaches to the problem of 
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power control in wireless networks have been investigated in 
|@), and surveyed in [6 |-|8|. In this paper, we are interested 
in the other scenario where the users are competing against one 
another and aiming to maximize their own information rates. 
This competitive rate-maximization problem can be modelled 
as a strategic noncooperative game. The Nash equilibrium |9| 
of this game can be achieved via a distributed waterfilling 
algorithm where each user performs waterfilling by consider- 
ing the multi-user interference as an additive coloured noise. 
However, most of the current results on rate-maximization 
and waterfilling algorithms assume the availability of perfect 
information which is a strong requirement and cannot be met 
by practical systems. 

This paper addresses the following fundamental questions: 
i) How can the users independently allocate power if the 
channel state information (CSI) they have is imperfect? How 
can we formulate a rate maximization game under channel 
uncertainty and what is the nature of the equilibrium of this 
game? ii) What are the existence and uniqueness properties of 
such an equilibrium? How can such a solution be computed 
by a distributed algorithm and what are the conditions for 
asymptotic convergence of such an algorithm? iii) How are 
these conditions affected by the channel uncertainty? iv) What 
is the effect of uncertainty on the sum-rate and the price 
of anarchy of such a system? In answering these questions, 
we can gain further insight into the behaviour of waterfilling 
algorithms and methods to improve sum-rate in general. 

A. Summary of Main Results 

The main contributions of this paper are three-fold. First, we 
provide a game-theoretic solution for the problem of competi- 
tive rate-maximization in the presence of channel uncertainty. 
Secondly, we analyse the efficiency of the equilibrium as 
a function of the channel uncertainty bound and prove that 
the robust waterfilling solution proposed achieves a higher 
sum-rate with increasing uncertainty under certain conditions. 
Finally, we verify these results via simulations. We show that 
improved sum-rate under channel uncertainty entails a cost in 
terms of more stringent conditions for the uniqueness of the 
equilibrium and slower convergence of the iterative algorithm 
to the equilibrium. 

For systems with bounded channel uncertainty, we present 
a distribution-free robust formulation of the rate-maximization 
game based on an ellipsoid approximation of the uncertainty in 
the CSI. We present the robust-optimization equilibrium (RE) 
for this game. At the equilibrium, the users perform a modified 
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waterfilling operation where frequency-overlap among users 
is penalized. We derive sufficient conditions for the existence 
and uniqueness of the equilibrium and for the convergence 
of an asynchronous iterative waterfilling algorithm to the 
equilibrium of this game. 

In our work, we investigate the effect of channel uncertainty 
on the sum-rate of the system for the two-user case under two 
scenarios, viz., a two-frequency system and a system with 
large number of frequencies. For the two-frequency system, 
the equilibrium sum-rate improves and the price of anarchy 
decreases as the channel uncertainty increases under high 
interference. On the other hand, the behavior is reversed at 
low interference. Another important contribution of our work 
is to show that increasing channel uncertainty always drives 
the equilibrium closer to a frequency division multiple access 
(FDMA) solution for any set of channel coefficients for the 
system with asymptotic number of frequencies. This is because 
the users become more conservative about causing interference 
under increased uncertainty and this leads to better partitioning 
of the frequencies among the users. Under certain channel 
conditions, this also translates to an improvement in the sum- 
rate and a decrease in the price of anarchy of the system. Thus, 
we show an interesting phenomenon where increased channel 
uncertainty can lead to a more efficient equilibrium, and hence, 
a better sum-rate in certain multi-user communication systems. 

B. Related Work 

An iterative waterfilling algorithm for maximizing infor- 
mation rates in digital subscriber line systems iflOl is an 
early application of a game-theoretic approach to designing a 
decentralized algorithm for multi-user dynamic power control. 
This framework has been further analyzed and extended in 
irTTI - |fT5l . The inefficiency of the Nash equilibrium (which 
need not be Pareto optimal) has been addressed in [16| and 
lUTI and methods to improve the sum-rate of the system by 
using various pricing schemes and modified utility functions 
have been presented in ifTSl . JTHJ, fl9l . A centralized controller 
maximizing the sum-rate of the system leads to a non-convex 
optimization problem and has been shown to be strongly NP- 
hard in ||20l . The Pareto-optimality of the FDMA solution 
for this sum-rate maximization problem under certain channel 
conditions has been proved in fl2~Tl . 

Uncertainty in game theory and distributed optimization 
problems has only recently been investigated. The issue of 
bounded uncertainty in specific distributed optimization prob- 
lems in communication networks has been investigated in [22 1 
wherein techniques to define the uncertainty set such that they 
can be solved distributively by robust optimization techniques 
are presented. In [23], incomplete-information finite games 
have been modelled as a distribution-free robust game where 
the players use a robust optimization approach to counter 
bounded payoff uncertainty. This robust game model also 
introduced a distribution-free equilibrium concept called the 
robust-optimization equilibrium on which our approach is 
based. 

A brief look at a robust optimization approach for the 
rate-maximization game with uncertainty in the noise-plus- 
interference estimate has been presented in 11241 . where the 



authors present a numerically computed algorithm unlike the 
closed form results presented here. Such a numerical solution 
prevents further mathematical analysis of the equilibrium and 
its behaviour under different uncertainty bounds. Also, this 
uncertainty model is different from ours, where we assume 
the availability of CSI of the interfering channels and that 
these quantities have a bounded uncertainty. 

A similar problem of rate-maximization in the presence 
of uncertainty in the estimate of noise-plus-interference lev- 
els due to quantization in the feedback channel has been 
considered in [25|. This problem has been solved using a 
probabilistically constrained optimization approach and as in 
our work, also results in the waterfilling solution moving closer 
to an FDMA solution, with corresponding improvement in 
sum-rate. However, the effect of quantization on the conditions 
for existence and uniqueness of the equilibrium and conver- 
gence of the algorithm have not been considered. The results 
presented in ll25l are for a sequentially updated algorithm 
whereas our results allow asynchronous (and thus sequential 
or simultaneous) updates to the algorithm. Also, the power 
allocations computed by such a probabilistic optimization for- 
mulation do not guarantee that the information rates expected 
will be achieved for all channel realizations, unlike our worst- 
case optimization formulation. Furthermore, the relative error 
(and not just the absolute error due to quantization) in the 
interference estimate as defined in [25 ] is assumed to be 
bounded and drawn from a uniform distribution, which is 
inaccurate. In addition, this bound on the relative error can 
only be computed if the noise variance at the receivers is 
assumed to be known (which is not the case). The bounds 
computed in such a fashion are very loose and will degrade 
system performance. The other assumption that this relative 
error bound is in the range [0, 1) means that the absolute 
quantization error has to be less than the noise variance at the 
receivers, which restricts the applicability of the approach. Our 
problem formulation has no such limitation on the uncertainty 
bound based on the noise variances in the system. 

Robust rate-maximization for a cognitive radio scenario 
with uncertainty in the channel to the primary user has 
been presented in [26 1 . This leads to a noncooperative game 
formulation without any uncertainties in the payoff functions 
of the game (unlike in our case) with robust interference limits 
acting as a constraint on the admissible set of strategies. This 
game is then solved by numerical optimization as there is no 
closed form solution. 

Initial results of our work, involving the problem formu- 
lation of the robust rate-maximization game and the results 
for sufficient conditions for existence and uniqueness of the 
equilibrium, along with the conditions for asymptotic conver- 
gence of the waterfilling algorithm to the equilibrium, have 
been presented in 0"). Initial results on the analysis of the 
efficiency of the equilibrium was presented in [2|. We now 
extend this work to include detailed proofs and further analysis 
of the efficiency of the equilibrium as a function of the channel 
uncertainty bound. 
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C. Paper Outline 

This paper is organized as follows: Section [TT] describes 
the system model and provides the necessary preliminaries. 
Section [Til] formulates the robust game model for the rate- 
maximization game for the single-input single-output (SISO) 
frequency-selective Gaussian interference channel. Section [IV] 
presents the robust-optimization equilibrium for this game 
and conditions for its existence and uniqueness, along with 
the conditions for asymptotic convergence of the waterfilling 
algorithm to the equilibrium. Section [V] presents the analysis 
on the effect of uncertainty on the sum-rate of the system for 
the two-user scenario. Section [VI] presents simulation results 
and Section IVIII presents the conclusions from this work and 
possible future research directions. 

II. System Model and Preliminaries 

Notations used: Vectors and matrices are denoted by lower- 
case and uppercase boldface letters respectively. The operators 
(•) T , (-) _T , E{-} and || • H2 are respectively transpose, trans- 
pose of matrix inverse, statistical expectation and Euclidean 
norm operators. The diagonal matrix with the arguments as 
diagonal elements is denoted by Diag(-). The quantity [A]jj 
refers to the (z, j)-th element of A. R™ x ™ is the set of 
m x n matrices with real non-negative elements. The spectral 
radius (largest absolute eigenvalue) of matrix A is denoted 
by p(A). The operation [x] b a is defined as [x] b a = a if x < 
a; x if a < x < b; b if x > b and (x) + = max(0,a;). 
The matrix projection onto the convex set i? is denoted by 
[x].g = argmin z6 ^ ||z — x||2. The term w > indicates that 
all elements of w are positive and X >- indicates that the 
matrix X is positive definite. 

We consider a system similar to the one in ifPfl . which is a 
SISO frequency-selective Gaussian interference channel with 
N frequencies, composed of Q SISO links. f2 = {1, . . . , Q} is 
the set of the Q players (i.e. SISO links). The quantity H rq (k) 
denotes the frequency response of the fc-th frequency bin of 
the channel between source r and destination q. The variance 
of the zero-mean circularly symmetric complex Gaussian 
noise at receiver q in the frequency bin k is denoted by 
i7g(fc). The channel is assumed to be quasi-stationary for the 
duration of the transmission. Let er 2 (fc) = a q (k)/\H qq (k)\ 2 
and the total transmit power of user q be P q . Let the vector 
s q = [s q (l) s q (2) . . . s q (N)} be the N symbols transmitted by 
user q on the N frequency bins and p q (k) = E{|s g (fc)| 2 } be 
the power allocated to the fc-th frequency bin by user q and 
Pq — [p<j(1) ■ • -Pq(N)} be the power allocation vector. 

The power allocation of each user q has two constraints: 

• Maximum total transmit power for each user 



N 



E{l|s g ||^}=^P 9 (fc)<P g 



(1) 



fc=i 



for q = 1, . . . , Q, where P q is power in units of energy 
per transmitted symbol. 
Spectral mask constraints 

E{| Sg (fc)| 2 }= Pg (fc)<p™-(fc), (2) 



for k = 1, . . . , N and q = 1, . . . , Q, where p q nax (k) is 
the maximum power that is allowed to be allocated by 
user q for the frequency bin k. 
Each receiver estimates the channel between itself and all 
the transmitters, which is private information. The power allo- 
cation vectors are public information, i.e. known to all users. 
Each receiver computes the optimal power allocation across 
the frequency bins for its own link and transmits it back to the 
corresponding transmitter in a low bit-rate error-free feedback 
channel. Note that this leads to sharing of more information 
compared to other works in literature such as [14]. The channel 
state information estimated by each receiver is assumed to 
have a bounded uncertainty of unknown distribution. Ellipsoid 
is often used to approximate unknown and potentially com- 
plicated uncertainty sets |27l . The ellipsoidal approximation 
has the advantage of parametrically modelling complicated 
data sets and thus provides a convenient input parameter 
to algorithms. Further, in certain cases, there are statistical 
reasons leading to ellipsoidal uncertainty sets and also results 
in optimization problems with convenient analytical structures 
EHl, (29]. 

We consider that at each frequency, the uncertainty in the 
channel state information of each user is deterministic ally 
modelled under an ellipsoidal approximatiorQ 

T q = <F rq {k) + AF rq>k : lA^fcl 2 < e 2 

r#9 (3) 

V* = 1 jvj 

where e q > V q G 51 is the uncertainty bound and 



F rq {k) 



\H r q{k)\ 



(4) 



r9VV \H qq {k)\ 2 ' 

with F rq (k) being the nominal value. We can consider un- 
certainty in F rq {k) instead of H rq (k) because a bounded 
uncertainty in Fr q (k) and H rq (k) are equivalent, but with 
different bounds |3 

The information rate of user q can be written as [ 30 1 

N 



i? 9 = J>g 1 



k=l 



Pq(k) 



, (5) 



where a q (k) = a q (k)/\H qq (k)\ 2 . The two popular measures 
of "inefficiency" of the equilibria of a game are the price of 
anarchy and the price of stability. The price of anarchy is 
defined as the ratio between the objective function value at 
the socially optimal solution and the worst objective function 
value at any equilibrium of the game ll3D . The price of 
stability is defined as the ratio between the objective function 
value at the socially optimal solution and the best objective 
function value at any equilibrium of the game [ 3 1 1 . We 
consider the sum-rate of the system as the social objective 

'More specifically, we have a spherical approximation in (5). 

2 The model considered here has some redundancy in the uncertainty for 
the case when F rq (k) = which leads to including F rq (k) + A_F r9 ^ < 
in the model which can never happen in practice. However, this does not 
affect the solution in our method due to the nature of the max-min problem 
formulation in jilt which leads to selection of positive values of AF rq ^. 
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function. The sum-rate of the system is given by 

Q 

s = EX (6) 

9=1 

In our case, the price of stability and anarchy are the same 
as we prove the sufficient conditions for the existence of a 
unique equilibrium in Theorem [T] Thus, the price of anarchy, 
PoA, is the ratio of the sum-rate of the system at the social 
optimal solution, S*, and the sum-rate of the system at the 
robust-optimization equilibrium, S roh , i.e., 

(7, 



PoA = 



5-rob 



Note that a lower price of anarchy indicates that the equilib- 
rium is more efficient. 



III. Robust Rate-Maximization Game Formulation 
A. Nominal Game - No Channel Uncertainty 

The problem of power allocation across the frequency bins 
is cast as a strategic noncooperative game with the SISO links 
as players and their information rates as pay-off functions 1 14 1. 
Mathematically, the nominal game <g nom can be written as, 

Vq G Q, 



max 

Pg 



N 

E 

fe=i 



log 1 



Pq(k) 



<J 2 q (k)+Y,F rq {k) Pr {k)' (8 ) 

r=£q 



S. t. p q G V, 



<f 



where Q = {1, . . . , Q} is the set of the Q players (i.e. the 
SISO links) and V q is the set of admissible strategies of user 
q, which is defined as 



V a 



Pg e 



V<P q {k)<p™ ax {k), 



N 

E 

fe=i 



Pq{k) = Pq, k = l,. 



.N 



(9) 



The inequality constraint in (Q~|i is replaced with the equality 
constraint in (0 as, at the optimum of each problem in (O, 
the constraint must be satisfied with equality Ifl4l . To avoid 
the trivial solution p q (k) = p q nax (k) Vfc, it is assumed that 

12k=iPq lax > Pq- Further, the players can be limited to pure 
strategies instead of mixed strategies, as shown in ifTTl . 



B. Robust Game - With Channel Uncertainty 



According to the robust game model 11231 . each player 
formulates a best response as the solution of a robust (worst- 
case) optimization problem for the uncertainty in the payoff 
function (information rate), given the other players' strategies. 
If all the players know that everyone else is using the robust 
optimization approach to the payoff uncertainty, they would 
then be able to mutually predict each other's behaviour. The 
robust game gf rob where each player q formulates a worst-case 



robust optimization problem can be written as, V q G f2 



max mm 



S.t. p„ G V q , 



E lo § 1 



Pq(k) 



° 2 q (k) + E F rq (k)p r (k) 

T=£q 



(10) 

where T q is the uncertainty set which is modelled under 
ellipsoid approximation as shown in ([3J. This optimization 
problem using uncertainty sets can be equivalently written in 
a form represented by protection functions 11221 as, \lq G O, 



JV 



max mm 



in ^ log I 1 



fe=l 



Pq(k) 



Pq AF rq , k ^ ~\ ^(k) + J2(Fr,lk) + AF rq . k )p r (k) 



r^q 



s.t. Y,\AF rq , k \ 2 < e 2 q , p q eV, 

r=jtq 

From the Cauchy-Schwarz inequality 



(ID 



we get 



Y, AF r q ,kPr(k)< ^|AF r9 , fe | 2 ^| Pr (fc)| 

r^q r^q r^q 



(12) 



Using ( fl2l i. we get the robust game gf rob as, \/q G O, 

Pq (k) 

-Ug, I j. -i 

k=l 



N 



max 

Pg 



E log ( 1 



o- 2 q(k)+j:F rq (k)p r (k)+eq jEpm 



r¥=q 



r^q 



S.t. p q G V q . 



(13) 



Now that we have defined the problem for robust rate- 
maximization under bounded channel uncertainty, we present 
the solution to the optimization problem in ( fT3l for a single- 
user in the following section. 



C. Robust Waterfilling Solution 

The closed-form solution to the robust optimization problem 
in ( fT3l for any particular user q is given by the following 
theorem: 

Lemma 1. Given the set of power allocations of other users 
P q = {Pii ■ ■ ■ iPg-liPg+lj ■ • ■ >Pq}> the solution to 



N 



max \^ log ( 1 + 
p 9 ^— ' » 



fc=i 



Pq(k) 



° 2 q{k)+Y.Frq{k)p r {k)+e q lY J P 2 r{k) 



r^q 



r^q 



S.t. p q G V q . 

is given by the waterfilling solution 

p q = RWF q (p_ g ), 



(14) 



(15) 



where the waterfilling operator RWF q (-) is defined as 
[RWF q (p_ 9 )] fe 4 



.-ipr x ( fc ) 



fx q -a 2 q (k)-^F rg (k) Pr (k)-e q Y,Pr( k ) 

r=jiq V r^q 



(16) 
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Input: 

fi: Set of users in the system 
V q : Set of admissible strategies of user q 
T q : Set of time instants n when the power vector p q n ' of 
user q is updated 

T: Number of iterations for which the algorithm is run 
T^(n): Time of the most recent power allocation of user r 
available to user q at time n 
RWF q (-): Robust waterfilling operation in ( [ToT l 
Initialization: n = and p q ^ <— any p e V q , Vg S O 
for n = to T do 



end for 



RWF q 
p (n) 



(r'(n)) 



if n e T q , 
otherwise, 



Vg e ft. 



for k = 1,...,N, where p, q is chosen to satisfy the power Algorithm 1 - Robust Iterative Waterfilling Algorithm 
constraint Ylk=lPq(k) ~ Pq- 

Proof: This can be shown using the Karush-Kuhn-Tucker 
optimality conditions ||27l of this problem. ■ 

The robust waterfilling operation for each user is a dis- 
tributed worst-case optimization under bounded channel un- 
certainty. Compared with the original waterfilling operation in 
lfT4l under perfect CSI (i.e. e q = 0), we see that an additional 
term has appeared in ( TToT ) for e q > 0. This additional term can 
be interpreted as a penalty for allocating power to frequencies 
having a large product of uncertainty bound and norm of the 
powers of the other players currently transmitting in those 
frequencies. This is because the users assume the worst-case 
interference from other users and are thus conservative about 
allocating power to such channels where there is a strong 
presence of other users. 

Having derived the robust waterfilling solution for a single 
user in the presence of channel uncertainty, we consider 
whether a stable equilibrium for the system exists and if so, 
what its properties are and how it can be computed in the 
multi-user scenario in the following section. 

IV. Robust-Optimization Equilibrium 

The solution to the game £f rob is the robust-optimization 
equilibrium (RE). At any robust-optimization equilibrium of 
this game, the optimum action profile of the players {p*} ge n 
must satisfy the following set of simultaneous waterfilling 
equations: Vq G ft, 



and 



[S' 



max F ra (k), 
kev„nv r q 

0, 



if r ^ q, 



otherwise, 



(20) 



p* = RWF q (pl, . . . , pj_ x , p* +1 , . . . ,p£) = RWF q (p* ,). 

(17) 

It can easily be verified that the RE reduces to the Nash 
equilibrium of the system [14| when there is no uncertainty in 
the system. In Section [V] we analyse the global efficiency of 
the RE and show that the RE has a higher efficiency than the 
Nash equilibrium due to a penalty for interference which en- 
courages better partitioning of the frequency space among the 
users. The robust asynchronous iterative waterfilling algorithm 
for computing the RE of game £f rob in a distributed fashion 
is described in Algorithm Q] 

A. Analysis of the RE of Game Sf rob 

Let M = {1, . . . , N} be the set of frequency bins. Let V° 
denote the set of frequency bins that user q would never use 
as the best response to any set of strategies adopted by the 
other users, 



2>° = {kE {!,... ,N) 



[RWF q (p_ 9 )] fe = yp- q eP-q}, 



(18) 



where V- q = V\ x • • • x V q -i x Vq+i x • • • x Vq. The non- 



negative matrices E and S r 



[El 



eg, 



qr 



0. 



are defined as 

if r ^ q, 
otherwise, 



(19) 



where V q is any subset of {1, . . . , N} such that N — T>° C 

■l\ {1 V}. 

The sufficient condition for existence and uniqueness of 
the RE of game ^ Iob and for the guaranteed convergence of 
Algorithm Q] is given by the following theorem: 

Theorem 1. Game {# rob has at least one equilibrium for 
any set of channel values and transmit powers of the users. 
Furthermore, the equilibrium is unique and the asynchronous 
iterative waterfilling algorithm described in Algorithm^ con- 
verges to the unique RE of game ^ Iob as the number of 
iterations for which the algorithm is run, T — > oo for any 
set of feasible initial conditions if 

p(S rnax ) < 1 -p(E), (21) 

where E and S are as defined in H9\ and ( 1201 l respectively. 

Proof: See Appendix lei ■ 
In the absence of uncertainty, i.e. when e q = Vq G f2, 
we can see that this condition reduces to condition (CI) in 
OH as expected. Since p(E) > 0, the condition on S max 
becomes more stringent as the uncertainty bound increases, 
i.e. the set of channel coefficients for which the existence of 
a unique equilibrium and the convergence of the algorithm is 
guaranteed shrinks as the uncertainty bound increases. Also, 
from Lemma[3] we can see that the modulus of the waterfilling 
contraction increases as uncertainty increases. This indicates 
that the convergence of the iterative waterfilling algorithm 
becomes slower as the uncertainty increases, as shown in 
simulation results. 

Corollary 1. When the uncertainties of all the Q users 
are equal (say e), the RE of the game Sf rob is unique and 
Algorithm Q] converges to the unique RE of game 
T — » oo for any set of feasible initial condition if 

p(S max ) < 1 - e(Q - 1) 



?rob 



(22) 



Proof: When the uncertainties of all Q users is e, we get 
p(E) = e(Q-l). ' ■ 

The above corollary explicitly shows how the uncertainty 
bound and the number of users in the system affect the 
existence of the equilibrium and the convergence to the 
equilibrium using an iterative waterfilling algorithm. For a 
fixed uncertainty bound, as the number of users in the system 
increases, there is a larger amount of uncertain information 
in the system. Hence, the probability that a given system for 
a fixed uncertainty bound will converge will decrease as the 
number of users in the system increases. Also, if e(Q— 1) > 1, 
we will not have a guaranteed unique equilibrium and algorith- 
mic convergence for non-zero uncertainty bounds regardless of 
the channel coefficients. This will help designers plan for the 
appropriate uncertainty bounds based on the planned number 
of users in the system. 



Tx 1 



[1,1] 



Rx 1 



[ma, a] 



Tx 2 



Rx 2 



[1,1] 



Fig. 1: Anti-symmetric system with Q = 2, N — 2, t\ = e 2 = 
e and the noise variances for both users in both frequencies 
is a 2 . The channel gains are: |i?n(l)| 2 = \Hn(2)\ 2 = 
\H 22 (1)\ 2 - \H 22 (2)\ 2 = l;|#i 2 (2)| 2 = |# 2 i(l)| 2 = a 
and |#i 2 (l)| 2 = |iJ 2 i(2)| 2 = ma with m > 1 and 
< a < 1. The power allocations for this system at the 
robust-optimization equilibrium are presented in 



V. Efficiency at the Equilibrium - Two-User Case 

In this section, we analyse the effect of uncertainty on the 
social output of the system. For the two user case, the worst- 
case interference in each frequency reduces to ( y F rq (k) + 
eq)p r {k) with q,r — 1,2 and q ^ r. This means that the 
robust waterfilling operation for the two user case (Q = 2) is 
simply the standard waterfilling solution with the worst-case 
channel coefficients. We restrict the analysis to the two-user 
case with identical noise variance <r 2 (fc) = a 2 V7c,<7 across 
all frequencies and identical uncertainty bounds e\ = e 2 = e 
and total power constraints J2k=i Pi(&) = J2k=i Pa(A) = Pt 
for both users. These results can be extended to the non- 
identical case along similar lines. In order to develop a clear 
understanding of the behaviour of the equilibrium, the sum- 
rate of the system is first analyzed for a system with two 
frequencies (N = 2) and then extended to systems with large 
(N —> 00) number of frequencies. 

A. Two Frequency Case (N = 2) 

Consider a two-frequency anti-symmetric system as shown 
in Figure Q] where the channel gains are |Hn(l)| 2 = 
\H n (2)\ 2 = \H 22 (l)\ 2 = \H 22 (2)\ 2 = l;| J ff 12 (2)| 2 = 
|# 2 i(l)| 2 = a and |#i 2 (l)| 2 = |if 2 i(2)| 2 = ma with m > 1 
and < a < 1. From ( TToT l the power allocations at the robust- 
optimization equilibrium of this system are, 

Pl (l) = [Vi - a 2 -( a + e)p 2 {l)]+, 
p 1 (2) = [ f i 1 -a 2 -(ma + e)p 2 (2)} + , 
P 2 (l) = L"2 - er 2 - (ma + e)pi(l)] + , 
p 2 (2) = [M 2 -^ 2 -(a + e)pi(2)] + , 

and the total power constraint for each user is pi(l) +pi(2) = 
p 2 (l) + p 2 (2) = 1. Let pi(l) = p, hence, by symmetry, 
Pi (2) = p 2 (l) = 1 -p, p 2 {2) = p and m = /i 2 = fi. 

The following theorem presents the effect of uncertainty on 
the sum-rate and price of anarchy of the system for the high 
interference and low interference cases. 

Theorem 2. For the two-user two-frequency anti-symmetric 
system described above, we have the following results: 



• High interference: When a 2 ^ a(l — p), the sum-rate 
increases and the price of anarchy decreases as the 
channel uncertainty bound increases. 

• Low interference: When a 2 3> map, the sum-rate de- 
creases and the price of anarchy increases as the channel 
uncertainty bound increases. 

Proof: See Appendix ITJ1 ■ 
We can see from this result that the RE behaves in opposite 
ways when there is high interference and when there is low 
interference in the system. This suggests that there might be 
a certain level of interference where the sum-rate and price 
of anarchy do not change with change in uncertainty. This is 
given by the following proposition: 

Proposition 1. At the level of interference a — a , 

where 

— £((<-' +i > 2+ !?) 4 -"'- 1 )' <24 » 

the sum-rate and the price of anarchy are independent of the 
level of uncertainty or the power allocation used. Furthermore, 
at this value of interference, the price of anarchy is equal to 
unity. 

Proof: See Appendix |E] ■ 
We can see that even for such a simple system, the global 
behaviour of the robust-optimization equilibrium appears to 
be quite complex. This indicates that the global properties 
of the robust-optimization equilibrium for larger systems is 
quite strongly dependent on the level of interference in the 
system, which is seen in the following results with asymptotic 
number of frequencies. However, the underlying nature of the 
algorithm for the two-frequency case is that the RE moves 
towards the FDMA solution as the uncertainty bound increases 
(from (|63]0. 

B. Large Number of Frequencies (N — > 00) 
We consider J(fc), defined as 

J(k) 4 - Pl (k)p 2 (k), (25) 
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as a measure of the extent of partitioning of the frequency k. 
It is minimum (J(k) = — 1) when both the users allocate all 
their total power to the same frequency k and is maximum 
(J(k) — 0) when at most one user is occupying the frequency 
k. Note that J(k) = Vfc G {1, . . . , iV} when the users adopt 
an FDMA scheme. 

The following lemma describes the effect of the uncertainty 
bound on the extent of partitioning of the system: 

Lemma 2. When the number of frequencies, N — > oo, the 
extent of partitioning in every frequency is non-decreasing as 
the uncertainty bound of the system increases for any set of 
channel values, i.e., 

d 

— J(k)>0 Vfc G {1, . . . , N} when N -> oo, (26) 
oe 

with equality for frequencies where J(k) — 0, where J(k) is 
defined in d25l > 

Proof: See Appendix [F] ■ 
The above lemma suggests that the robust-optimization 
equilibrium moves towards greater frequency-space partition- 
ing as the uncertainty bound increases when there is a large 
number of frequencies in the system. In other words, the RE is 
moving closer to an FDMA solution under increased channel 
uncertainty. When the FDMA solution is globally optimal, 
this will lead to an improvement in the performance of the 
equilibrium. This is stated in the following theorem: 

Theorem 3. As the number of frequencies, N — >■ oo, the sum- 
rate (price of anarchy) at the robust- optimization equilibrium 
of the system is non-decreasing (non-increasing) as the uncer- 
tainty bound increases if, V k G {1, . . . , N}, 

[F 21 {k) - e)(F 12 (k) - e) > ~. (27) 

Proof: Using [21 1 Corollary 3.1], we find that the sum of 
the rates of the two users in the frequency k is quasi-convex 
only if F2i(k)Fi2(k) > 1/4. Let C be the minimum number 
of frequencies occupied by any user. When there are only 
two users and a large number of frequencies, C > 1. If the 
condition F2i(k)Fi2(k) > j(l + tj^tj) 2 i s satisfied for some 
C > 2 for all frequencies k G {1,...,N} (thus satisfying 
F2i(k)Fm(k) > 1/4), then the Pareto optimal solution is 
FDMA ||2T1 Theorem 3.3]. This needs to be satisfied for the 
worst-case channel coefficients which leads to $T7\ . Thus, the 
solution moving closer to FDMA will improve the sum-rate 
of the system. From Lemma [2] the robust equilibrium moves 
closer to FDMA as uncertainty increases and thus will result 
in an improvement in sum-rate. 

The Pareto optimal solution under this condition (which 
is FDMA) is constant under varying uncertainty bounds as 
such an uncertainty in the interference coefficients Fi2(fe) and 
i<2i(fe) does not affect the FDMA solution where there is no 
interference. Thus, an increase in sum-rate will result in an 
decrease in price of anarchy. ■ 

Remark 1. For the special case of frequency-flat systems, at 
the equilibrium, all users have equal power allocation to all 
frequencies and this is not dependent on the uncertainty in 
the CSI. This leads to no change in the extent of partitioning 



and thus sum-rate and price of anarchy are not affected by 
uncertainty. 

Remark 2. The results of this section are not just limited to the 
robust-optimization equilibrium for the system presented here. 
When the uncertainty e = 0, the framework presented here can 
be used to analyse the behaviour of the Nash equilibrium of the 
iterative waterfilling algorithm as a function of the interference 
coefficients. 

Remark 3. The modified waterfilling operation in ( TToT ) can also 
be used as a pricing mechanism to achieve improved sum-rate 
performance in a system with no uncertainty where e is a 
design parameter, with all the analytical results presented here 
still being valid. 

VI. Simulation Results 

In this section, we present some simulation results to study 
the impact of channel uncertainty on the RE by comparing 
it with the ideal scenario of NE under perfect CSI. Figure |2] 
shows the simulation results for the two user and two fre- 
quency scenario and Figures l3Tfl0l show the results for a more 
general system under different settings. 

In Figure [2a] we can see the sum-rate at high interference 
as a function of interference and uncertainty in the system 
shown in Figure [T] The flat region corresponds to the sum- 
rate at Pareto optimal solution (FDMA) and the edge of 
the surface corresponds to the sufficient condition in $2~T[ . It 
can be seen that the Nash equilibrium (uncertainty, e = 0) 
moves closer to the Pareto optimal solution as the interference 
increases. It is also evident that the sum-rate increases for a 
fixed interference as uncertainty increases, as expected from 
Theorem [2] In Figure |2bJ we can see the sum-rate at low 
interference as a function of interference and uncertainty. 
As expected from Theorem [2] the sum-rate decreases as the 
uncertainty increases. 

In Figures I3TH01 the behaviour of the equilibrium under 
varying uncertainty bounds is investigated through numerical 
simulations. The simulations are computed for a system with 
Q users and N frequencies averaged over 5000 channel 
realizations. The channel gains are H rq (k) ~ Nc(0,l) for 
r ^ q and H qq (k) ~ Nc(0, 2.25). The channel uncertainty 
model has nominal value F rq (k) — F^ q (k)(\ + e rq (k)) with 
e rq (k) ~ U(— |, |), S < 1. The specific parameter values 
used for the simulations are provided above each figure. Note 
that the zero uncertainty solution corresponds to the Nash 
equilibrium and the nominal solution is the solution resulting 
from using the erroneous channel values in the traditional rate- 
maximization game & s without accounting for its uncertainty. 
The effect of uncertainty, number of users and number of 
frequencies on the average sum-rate of the system, the average 
number of frequencies occupied by each user and the average 
number of iterations for convergence are examined. In these 
figures, the Nash equilibrium point is when the uncertainty is 
zero. 

In Figure [3] it can be observed that the sum-rate at the 
Nash equilibrium under perfect CSI is less than the sum-rate 
at the robust-optimization equilibrium under imperfect CSI 
and that the gap in performance increases as the uncertainty 5 
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Fig. 2: Simulation results for the anti-symmetric system in 
Figure Q] Note that the zero uncertainty corresponds to the 
Nash equilibrium 
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Fig. 8: Average number of iterations vs. uncertainty, S. 
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Fig. 9: Average number of iterations vs. number of users, Q. 

increases. Under imperfect CSI, the power allocation using the 
robust-optimization equilibrium^ in ( fToT ) and dl7l > has higher 
sum-rate as uncertainty increases, because the users are more 
cautious about using frequencies with significant interference, 
thus reducing the total amount of interference in the system. 
The Pareto optimal solution in this scenario is some form of 
FDMA (where the specific channels allocated to a particular 
user will depend on the actual channel gains), and does not 
depend on the channel uncertainty bound for a given system 
realization. Thus, when the sum-rate of the system increases 
with rise in uncertainty, we can expect that the price of anarchy 
will decrease. 

In Figure |4j it can be observed that the sum-rate of the 
system under the robust solution reduces when the number of 
users increases. This is because having a greater number of 
users results in higher interference for all users and this effect 
is strong enough to counter user diversity which would have 
resulted in higher sum-rates if the users were on an FDMA 
scheme. In Figure[5] it can be observed that the sum-rate of the 
system improves with increase in number of frequencies and 
also that the robust solution continues to perform better than 
the nominal solution even when the number of frequencies 
increases. 

In Figure [6] it can be seen that the robust solution results 
in a lower average number of channels per user as the 
uncertainty, S increases. Also, the total number of channels 

3 {16} and (T7J are in terms of absolute uncertainty e while the simulations 
use relative uncertainty <5. They are equivalent to one another. 
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Fig. 10: Average number of iterations vs. number of frequen- 
cies, N. 



each user occupies at the robust-optimization equilibrium is 
lesser than at the nominal solution, regardless of the number 
of users, as can be seen in Figure [7] This implies that the 
users are employing a smaller number of frequencies, which 
demonstrates the improved partitioning of the frequency space 
among the users in order to reduce interference. Hence, this 
leads to the higher sum-rates as observed in Figure [3] 

In Figures [8] and [9] it can be observed that the average 
number of iterations for convergence increases as the uncer- 
tainty S, and the number of users Q increase respectively. 
This is as expected from Lemma [3] and Corollary Q] as the 
modulus of the block-contraction in ( f39l > increases as the 
uncertainty increases. This indicates that the step size of 
each iteration reduces as uncertainty increases, leading to 
slower convergence. In Figure \TU[ it can be observed that 
the average number of iterations increases with the number 
of frequencies. This is due to the fact that when there are 
more frequencies, and as such there is a greater probability 
of a channel realization being drawn from the tail of the 
Gaussian distribution used to generate them, which results in a 
smaller modulus of the block-contraction in d39b on average, a 
greater number of iterations are required in order to converge 
to the equilibrium. Thus, the trade-off for robust solutions with 
higher sum-rates is in a higher number of iterations before 
convergence. 

VII. Conclusions 

In this paper, we have presented a novel approach for rate- 
maximization games under bounded channel state information 
uncertainty. We have introduced a distribution-free robust for- 
mulation for the rate-maximization game. The solution to this 
game has been shown to be a modified waterfilling operation. 
The robust-optimization equilibrium (RE) for this game has 
been presented and sufficient conditions for its existence, 
uniqueness and asymptotic convergence of the algorithm to 
the RE have been derived. For the two-user case, the effect 
of uncertainty on the social output of the system has been 
analyzed. We analytically prove that the RE moves towards 
an FDMA solution as the uncertainty bound increases, when 
the number of frequencies in the system becomes asymptot- 
ically large. Thus, an interesting effect of improvement in 
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sum-rate as the uncertainty bound increases is observed. In 
summary, for systems with significant interference, bounded 
channel uncertainty leads to an improved sum-rate but at the 
cost of greater number of iterations. This framework can be 
extended to MIMO rate-maximization games, cognitive radio 
with various interference constraints and other noncooperative 
games. 



Appendix 

A. Robust Waterfilling as a Projection Operation 

Let Q q (k) represent the denominator terms in (TBI , which 
is the worst-case noise+interference 



4 a 2 (fc) +Y / F rq (k) Pr (k) + e q (28) 

r=£q V r^q 



It has been shown in [33] that the waterfilling operation 
can be interpreted as the Euclidean projection of a vector 
onto a simplex. Using this framework, the robust waterfilling 
operator in ([Tol l can be expressed as the Euclidean projection 
of the vector $ g = [$ g (l), . . • , $ q (N)] T onto the simplex V q 
defined in (0: 



RWF q (p_,) = [-*,] Pi 
which can be conveniently written as 

RWF q (p_,) = 



where 



rq 



(7 q ^ F r gP r ^q^q 

r^q 



4 [**(!),,.., * 2 q (N)] T , 
4 Diag(F r ,(l),...,F rff (2V) 



(29) 
(30) 

(31) 
(32) 

■ (33) 



B. Contraction Property of the Waterfilling Projection 

Given the waterfilling mapping RWF(-) defined as 

RWF(p) = (RWF q (p_,)), en ■ V^V, (34) 

where V = Vi X •■• x Vq, with V q and RWF q (p_ 9 ) 
respectively defined in (0 and (13 Oi l, the block-maximum norm 
is defined as ll34l 



RWF(p)| 



2, block 



max • 

gen 



RWF q (p 9 )| 



(35) 



where w = [w\, . . . , wq] t > is any positive weight vector. 
The vector weighted maximum norm is given by ll32l 



||x||- vec ^max^, w>0, xel« 
qen w q 

The matrix weighted maximum norm is given by 11321 
1 



|[A] gr .|w r , 



A G 



(36) 



(37) 



The mapping RWF(-) is said to be a block-contractiorQ 
with modulus a with respect to the norm 1 1 • 1 1 bl k if there 
exists a G [0, 1) such that, Vp (1) , p (2) G V, 

1 1 RWF(p( 1 ') -RWF(p( 2 ))| |" Wock < a| |p (1) -p (2) | |" block , 

(38) 

where pW = (p q l \ . . . , p W) for i = 1,2. 

The contraction property of the waterfilling mapping is 
given by the following lemma: 

Lemma 3. Given w = [wi,...,wq] t > 0, the mapping 
RWF(-) defined in (1341 l satisfies 

||RWF(p«) -RWF(p( 2 ))|| w .. , < 

w X I In' 
oo,mat M r 1 



1 2, block 
,(1) 



n (2)|| w 
v II 2, block' 



(39) 



Vp^\p (2) G V, where E and S as defined in (O and (l20b 
respectively. Furthermore, if 



El 



- < 1, (40) 



for some w > 0, then the mapping RWF(-) is a block 
contraction with modulus a = IIS""™ + Ell™ mn t. 



Proof: This proof is structured similar to that of 
Proposition 2], However the additional term e q { q in (l30l l 
necessitates a separate treatment here. 

(i) 

For each q G £1 and i = 1,2, given iq ' = 

Er^lW 5 . • ■ • . ^ let P-,(fc) W 4 

pi(k)M, . . . .Pj-iCAjW^^^AjW, . . . ,p Q (fc)W] and 
Af „A|| f W_ f (2)r 



Then. 



Af «= £L/£tf(*) (1 



< 




/^p2 (fc) (2) 



|P-,W (2) | 



fc— 1 r^q 



^l| P _ 9 (fc)w- P _ 9 (fc)( 2 ) 
fe=i 

AT 

EE(^( fc ) (1) +^( fc ) 

-2p r (fc)( 1 )p r (fc)< 2 ) 



(41) 



(42) 



(43) 



\(2) 



(44) 
(45) 



ElK-p? 3 

where (|43l follows from 11321 Lemma 5.1.2]. Now, define for 
each q € Q, 

RWF q (p-l) - RWF q (p^ q 



6RWF„ 



p(l) 



p (2) 



(46) 



4 The mapping T is called a block-contraction with modulus a 6 [0, 1) if 
it is a contraction in the block-maximum norm with modulus a 1341 . 
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Then, using ( f30l > in d46j l, eRWF q can be written as 



6RWF„ — 



-a -Vf d (1) -e f (1) 



r^q 



-[-^-E F ^ 2) - e ^ (2 



rjtq 



< 



£F rg pW+ £g f« 

rjtq 



E F ^ 2) 



e f ( 2 ) 
t q l q 



< 



e,K (1) -ff 



f (1) _ f (2) 
9 9 



(47) 

2 



(48) 



(49) 



(50) 



< E*% p^-pf 2) 



<E( m f ^,(*))||p^-p^ 



Elk^-p^ 



(51) 



E([S m0X ]r,+ 



(52) 

(53) 
(54) 



r^q 



V Pg^pf 5 £ and V g £ O, where: © follows from 
the nonexpansive property of the waterfilling projection [33 
Lemma 3]; d50b follows from the triangle inequality [32|; (l5TT i 
follows from (145k (l52l and (l53l follow from the definitions of 
F rq and e g respectively from (|3T1 > and (|46T >; and (l54l i follows 
from the definition of S max in (|20b and Jensen's inequality 

E3. 

The set of inequalities in d54"l > can be written in vector form 



as 



< e RW F < (S ? 



E)e 



(55) 



where E is defined in JT9l and the vectors enwF an d e 
are defined as e^wF — [grwFi > • - • > grwFq] , and e = 

[ei , eo\ T . Using the vector and matrix weighted maximum 

norms from (|3oT > and (|37| | respectively, d55l can be written as 



1 l w 


1 l e l 1 


1 1 oo,mat 





e RWF <US max +E)e 

I I vvv I I oo,vec — I I v 

< \ \S max +E 
V w > 0. Using the block-maximum norm (l35l l. we get 

||RWF (p (1) ) -RWF (p (2) 

<||S maa: + E|r 



(56) 



2, block I I ^RWF | | qq jVCC 



2, block 



(57) 



Vp(2) jP ( 2 ) e V, with E and S as defined in ([19) and (|20]l 
respectively. It is clear that RWF(-) is a block contraction 



when ||S— +E||- mat <l. 



C. Proof of Theorem Q] 

From (35 j, every concave gam^l has at least one equilib- 
rium. For the game $f rob ; 

1) The set of feasible strategy profiles, V q of each player 
q is compact and convex. 

2) The payoff function of each player q in ( fT3l ) is contin- 
uous in p G V and concave in p q G V q . 

Thus, the game 5f rob has at least one robust equilibrium. 

From Lemma [3j the waterfilling mapping RWF(-) is a 
block-contraction if (l40l > is satisfied for some w > 0. Thus, 
the RE of game Sf rob is unique (using [36, Theorem 1]). Since 
S ma;E +E is a nonnegative matrix, there exists a positive vector 
w such that 



E| loo, mat < 1 



(58) 



Using 11341 Corollary 6.1] and the triangle inequality [32], this 
is satisfied when 



^max 1 1 w 

' I loo, mat 



IEI 



oo,mat 



< l => P {s r 



< l-p(E). 

(59) 

From Lemma [3] and d58T > the waterfilling mapping RWF(-) 
is a block-contraction. From [36 , Theorem 2], the robust itera- 
tive waterfilling algorithm described in Algorithm Q] converges 
to the unique RE of game & Ioh for any set of feasible initial 
conditions and any update schedule. 

D. Proof of Theorem [2] 

Consider the interior operating points of the robust waterfill- 
ing operator RWF q ( ) where it is linear. Eliminating /1 from 
(l23l . we get 

V = W , ) ~ °~ ' r > 0.5. (60) 

2(1 — (to + l)a/2 — e) 

The signal-to-interference-plus-noise ratio (SINR) for the 
two users in the two frequency bins is given by 

P 



SINRi(l) = SINR 2 (2) = 
SINRi(2) = SINR 2 (1) = 



(T 2 + a(l — p) 

I-P 
a 2 + map 



(61) 



and the sum-rate of the system at the RE is 

S Ioh 4 2 log (l + P -) + 2 log (l , , 

v a + a[l — p)y \ a + map 



l-p 



The gradient of p with respect to e is 

dp (m — l)a 



> 0. 



(62) 
(63) 



de 4(1 - (m+ l)a/2 - e) 2 

Thus, the RE moves towards the FDMA solution as the 
uncertainty bound increases. 

5 A game is said to be concave if the payoff functions are concave and the 
sets of admissible strategies are compact and convex. 
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Case 1: High interference scenario: In the high interference 
scenario, a 2 <C a(l — p). Let £ = p/a(l — p). Then, the 
SINR for the two users in the two frequency bins can be 
approximated as 

P 



SINRi(l) = SINR 2 (2) 
SINRi(2) = SINR 2 (1) 



a(l -p) 
1-p 1 



(64) 



map ma 2 t; 

The sum-rate of the system at high interference can be 
approximated as 



S«21og((l+£)(l + — L-)) 
= 21og(l + -i- 7 +£ ' 



(65) 



ma' 



ma 2 !;' 



Our aim is to analyse the behaviour of the sum-rate S as 
the uncertainty e increases. To this end, we show that the 
gradient of the sum-rate with respect to e is positive. As log(x) 
increases monotonically with x, we consider 

(i-p)y c 

(66) 



9e ma 2 ^' 



(1 



1 



2i-n 



mp z 



(67) 



and (1 - ^ m pl ) > since p > 0.5 and m > 1. Now, 

d£ _ 1 dp 
de a(l — p) 2 de 

From (|63]l, d66]l and (j67]i, we get |f > 0. Thus, we see 
that the sum-rate of the system increases as the uncertainty 
e increases. This also shows that the robust-optimization 
equilibrium achieves a higher sum-rate in the presence of 
channel uncertainty (e > 0) than the Nash equilibrium at zero 
uncertainty (e = 0). 

The social optimal solution for this system at high interfer- 
ence is frequency division multiplexing 1211 . In other words, 
the frequency space is fully partitioned at the social optimal 
solution. The sum-rate at the social optimal solution for the 
given system at high interference, S*, is given by 

S* = 21og(l + 4)- (68) 



The price of anarchy at high interference, PoA, is 



PoA = 



log(l 



1 

in a - 



(69) 



Since |f > 0, we have ff^A < . 



Case 2: Low interference scenario: In the low inter- 
ference scenario, i.e. when map <C cr 2 , the signal-to- 
interference+noise ratio SINR for the two users in the two 
frequency bins can be approximately written as 

P_ 

(70) 

SINRi(2) = SINR 2 (1) 



SINRi(l) = SINR 2 (2) 



1 



V 



The sum-rate of the system at low interference can be approx- 
imated as 

1 



S«21og((l 
= 21og(U 



P_ 

a 2 
1 

^2 



)(1 



P 



P-P 



a- 



Now, 



dS_ 

de 



= 1 



1 



p - p 2 \ ~ l (1 - 2p) dp 



< 0. 



(71) 



(72) 



a" a" J a 2 de 

At low interference, the system behaves similar to a parallel 
Gaussian channel system. The social optimal solution in this 
scenario is the classical waterfilling solution and leads to equal 
power allocation, i.e., = Pi (2) = J> 2 (1) = P2(2) = p = 

0.5. The sum-rate at the social optimal solution for the given 
system at low interference, S*, is given by 

S* = 41og(l + ^). (73) 

The price of anarchy at low interference, PoA, is 

41og(l 



PoA 



l 

2a 1 



21og(l + ^ 
log (1 + £ H 



v-v 



1 

TP- 



log(l 



Note that, at low interference, map <C 1. From 



(74) 



we 



get p ?» 0.5. Thus the PoA is close to unity. Since < 0, 



we have > o. 



E. Proof of Proposition [7J 

The gradient of the sum-rate 5 ro b with respect to e is 

8S rob dS roh dp 



di 



dp de 



(75) 



From d63l , we have ^ > 0. Now, 



ap 



d grob 2 <7 2 +a(l-p) ^ (a 2 +a(l-p)) 2 

dp 



(7 2 +a(l— p) 
1 



(1-P)r 



(76) 



a z +map (a 2 +map) 



1 



J 1 ~P 

(j 2 -\-map 



Setting Qp Tab = 0, we solve for a to get the following roots, 



0, 



= { ym + 1 ± [Am/ a 2 + (m + 1 

(m-l)p 2 +2p-l ■ 



(77) 



The positive root that is independent of e and p (which is a 
function of the uncertainty e, from (l60i i) is the required solu- 
tion where the sum-rate is constant regardless of uncertainty. 
Thus, the required interference value is given by 

Oo = |- (Um/a 2 + (m+ l) 2 )^ - m - l) . (78) 
Lm \ I 



Since the root a of 



dS IO 
dp 



is independent of p, different power 



allocation schemes (resulting in different values of p) will 
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result in the same sum-rate at a = a a . Thus, the price of 
anarchy at a — a a is unity. ■ 

F. Proof of Lemma [2] 

From ( [TBI , the power allocations for the two users at the 
robust-optimization equilibrium in the fc-th frequency are 



(79) 



where 



Pi(fc) = -a 2 - (F 21 + e)p 2 (k) 
P2(k) = (fi 2 - v 2 - (F12 + e)pi(k) 



Pi(k) + (F 2 i(k) + e)p 2 (k) - m - a 2 
(F 12 (k) +e) Pl {k) + p 2 (k) - fi 2 - a 2 

E fce p o! P2{k) + n 2 {p, 2 - a 2 ) 
Writing these in matrix form we get, 



0, fc€2? o ,(80) 
0, fcG2? o; (81) 



Pt, 
P T . 



A fel 













~o 2 







h 




P(koi) 




o 2 


I 2 


I 2 


D 








Pt 



(82) 
(83) 



(84) 



where 

P(*) = 
Let 



1 F 21 (k)+e 
F 12 (k)+e 1 

Pi(k) 
p 2 (k) 



D 



ni 
n 2 



A 


'Pt 


A 


'fix - a 2 ' 


, Pt = 


P t_ 


M = 


M2 - o- 2 _ 



(85) 





~A kl 







h 





















h 



(86) 



C±[I 2 ... I 2 ] andP^ [p(h)...p(k ol )]- 
so that we can write 



A 


B 




P 




0" 


C 


D 




A 




Pt. 



We can solve this system to get 



p 




A B 


-l 


o" 




w 


X 




o" 




x Pt " 






C D 




Pt. 




Y 


Z 




Pt. 




Zp t 



(87) 



(88) 



w 


X 


A^ 


A 


B 


Y 


Z 




C 


D 



(89) 



Using ||37l Fact 10.12.9] and differentiating ( 1881 with respect 
to e, we get 



d_ 

de 



A B 
C D 



4- A 



A B 
C D 



(90) 
(91) 



with = Ef=iP2(fc) = Pt- 

Let T>i , 2?2 and V Q i be the sets of frequencies exclusively 
used by user 1, user 2 and by both respectively and n% = \T>i\ 
and n 2 = \T> 2 \ be the number of frequencies exclusively used 
by user 1 and user 2 respectively at the equilibrium. Then, 
from ( |79l , we have p\(k) = fix — o\ and p 2 (k) = V k 6 
2?i and px(k) = and p 2 (k) = fi 2 - a\ V k 6 V 2 . The 
power remaining for allocation to the frequencies in 1> i = W 
{fci, . . . , k i} by user 1 and user 2 is (1— ni/ii) and (1 — n 2 fi 2 ) 
respectively. 

This separation of the frequency-space into exclusive-use 
and overlapped-use frequencies allows us to analyse the sys- 
tem without the nonlinear operation (•)+. Thus, we can write 
the power allocations at the fixed point in the overlapped-use 
frequency-space as a system of linear equations, 





"W(|A)Xp t ~ 
Y(!A)Xp t 

Due to the nature of the waterfilling function, n\ and n 2 
are non-decreasing piecewise-constant functions of e. The 
above derivative exists only in the regions where n\ and 
n 2 are constant. From ||37l Proposition 2.8.7] and using 
A" 1 = Diag(A^\ . . . , A"j), we get 



""Ar 1 - A.^ZA" 1 



X=[A- i 1 Z ... A-jZ] T . 
Y = [ZAr 1 



-A^ZA" 1 



Ar 1 - A^ZAr 1 



Z = ^ 



n 2 



za; 



Ex £ 



E A, ^ 



F 2 i(i)+e 



E 



i 

A" 



,(92) 

(93) 
(94) 

(95) 



where A, 4 det(A l ) = 1 - (F 2 i(i) + e) (F 12 {i) + e), 



«2 



E 



kev a 



and 



F 2 i(i)+e 



1 

A, 
Fl2(»)+e 
A, 



EfceX> o; Ai 

F 12 {i)+e 



F 2 i(i)+e 
A, 
1 

Ai 



Thus, from dM) and (j9T|i, we get 
P = X Pt = [A^Zp t 



A fe ~l Z Pt]' 



and 
dP 

~d7 



-W(^)X, 
^Aj^zVGA^Zpt -A^GA^Zp t 

i— ki 



EA^ZA^^VZpt -A^GA^Z Pt 



(96) 



(97) 



(98) 



(99) 



(100) 



where 



G 



9Ai 
de 



1 

1 



V* = 1,...,N. 



Therefore, for fc = k\, . . . , k a i, 

p(k) = [ Pl (k) P2 (k)] T = A-^pt, 



(101) 



(102) 
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p'(fc) = -p(k)=[ P > 1 (k) P ' 2 (k)\ 

koi 

= Aj^ 1 zX)Ar 1 GA 4 - 1 Zpt 



(103) 



i—k± 



-A^GA^Zpt (104) 



Consider the extent of partitioning J(k) at any frequency k. 
For frequencies where J(k) = 0, at least one of the users has 
zero power allocation and that will not change with change in 
uncertainty. Thus J(fc) = for k e ViL)V 2 . Also, in cases 
when n\ or n 2 change due to some frequency k dropping from 
the set V i, J(k) increases from some negative value to zero. 

Now consider the extent of partitioning for frequencies 
where both users have non-zero power allocation. Differen- 
tiating the extent of partitioning for frequency k G V a i with 
respect to e, we get 



d_ 

~de 



J(k) 



d 



de 



(-Pi(k)p 2 (k)) 



(105) 

= -p[(k)P2(k)+ Pl (k) P ' 2 (k) (106) 
= -p(fc) T Gp'(fc) (107) 

= -(A,T 1 Z Pt ) T G(-A,T 1 GA,T 1 Z Pt 

Ki 

+A^ l Zj2 A 7 lG K 1 'Z'Pt) (108) 



i—ki 



= Pt^A- GA-^I 

kol 

i=fei (109) 

Let q fe = GA^Zpt. Using G T = G 1 = G, GA^G = 
A^ T and GA fc G = AT we get 



0_ 

7k 



J{k) = q^A^ 1 - A^Z K l K T **)<ik (HO) 



i—ki 



Let M fe = E-=' fcl A^A^AJ and Q k = A^ 1 -A^ZM^. 
When n a i — o(N) (i.e., when limjv->oo = 0), we 
have the total number of frequencies, n\ + n 2 = O(N). 
Since A~ 1 A~ T A^ = 0(1) for each i and k, we have 
M fe = O{n oi ) and Z = 0(1/N). Thus, 



lim A^ZM 



N- 



lim Qfe + = A fc 

iV— J-OO 



Afe T >■ o 



(111) 



p2xl 



as its 



from Theorem Q] Thus we get x T QfeX > Vx £ 
symmetric part Qfe + QjT is positive definite l38l . Hence, we 
get §- e J(k) > when N oo, with equality when J(k) = 
0. ■ 
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